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ABSTRACT: Determining the mechanism of charge transport through native DNA 
remains a challenge as different factors such as measuring conditions, molecule 
conformations, and choice of technique can significantly affect the final results. In this 
contribution, we have used a new approach to measure current flowing through isolated 
double-stranded DNA molecules, using fullerene groups to anchor the DNA to a gold 
substrate. Measurements were performed at room temperature in an inert environment 
using a conductive AFM technique. It is shown that the z-stacked B-DNA structure is 
conserved on depositing the DNA. As a result, currents in the nanoampere range were 
obtained for voltages ranging between +1 V. These experimental results are supported by a 
theoretical model that suggests that a multistep hopping mechanism between delocalized 
domains is responsible for the long-range current flow through this specific type of DNA. 


E INTRODUCTION between neighboring base pairs. An additional complication in 


The elucidation of the charge-transport mechanism along DNA comparing different conductance experiments is the nature of 
the anchoring groups that plays a key role in the overall 


is of prime interest for the design of nanoelectronic i n > E DNA2Z2” Th F 
components !~* as well as for biological*® and medical electronic characteristics o . ese effects are 


applications. ~! Different experimental techniques such as partially due to their intrinsic electronic properties of the 
photoinduced charge transfer’* or single molecule conductance anchoring groups that dictate the energencs of the injection 
experiments’? can be used to probe the transfer of charges process. In addition, the way the anchoring groups connect the 
along a DNA molecule. In the latter the molecule is connected DNA molecule to the surface determines the magnitude of the 
to macroscopic electrodes that are used to establish a steady electronic coup ling to the electrode. These groups can for 
state current flowing through the molecular structure. This can example be intercalated in the SEquence’Or located at the end of 
be achieved either using break junction techniques or by the DNA molecule. Finally, the anchoring groups can also affect 
depositing the molecule on a substrate and probing its the integrity of the D NA structure, and the strength and 
conductance using the tip of an STM or an AFM. Using geometry of the anchoring groups can either stabilize or 
such conductance techniques, a variety of electrical behaviors destabilize the double-helical structure of the DNA and lead to 


have been reported for DNA molecules, ranging from an large variations in the transport through the molecule itself. 
2 
13—16 


f ; 17-22 : Numerous theoretical models have been developed to 
Waa 23,24 > DS 3 near near simulate the transport of charges along DNA sequences. *° 
conductor." These marked differences can partially be linked P 8 8 eqt i 

to the preservation of the double-helical B-DNA structure. It Early works suggested the existence of two distinct mecha- 
has been described that for DNA molecules deposited on a nisms—superexchange and incoherent hopping” —to 


explain the charge transport along DNA molecules. According 
to these reports, the former is responsible of the charge 
transport along short DNA sequences containing up to three 
base pairs™ while the latter dominates the propagation of 


surface showing apparent heights of about 2.4 nm in AFM 
measurements behave as conductors. This height roughly 
corresponds to the diameter of B-DNA in which there is a 
relatively ordered a-stacked pathway for charge transport. In 
contrast, DNA molecules with an apparent height of about 1.1 
nm behave like insulators, most likely due to the disruption of Received: January 12, 2017 
the B-DNA structure that diminishes the -stack interactions Published: January 17, 2017 
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charge along longer sequences.*””” Recently, the existence of 
intermediate charge transport mechanisms, mediated by either 
delocalized conduction channels****~** or flickering resonan- 
ces,"* has been proposed and to some extent observed 
experimentally.*'**~*° Despite their intrinsic differences, 
these recent models agree on the importance of thermal 
fluctuations within the DNA structure.*”** These fluctuations 
facilitate charge propagation by introducing quasi-resonance 
between delocalized domains formed by a few strongly 
interacting base pairs.“ 

Previous attempts to probe the electronic behavior of DNA 
locally while „preserving its native structure have proven to be 
challenging.” We propose here an alternative way to study 
charge transport through single DNA molecules in their native 
structure by using fullerenesî! as anchoring groups in 
combination with conductive-AFM (C-AFM) at ambient 
temperature. The use of fullerene groups is interesting due to 
their small size, electronic properties, high affinity for gold, and 
possibility to be used in single molecule experiments.” °' The 
strong binding of the fullerene to gold surface leads to strong 
electronic couplings that should improve the overall charge 
transfer mechanism. To our knowledge, this approach has not 
been reported before. Using this new technique, we observed 
long-range charge transport over more than 20 nm and with 
current intensities in the nanoampere range. The good 
transport properties of the DNA-—fullerene (DNA-—C61) 
complexes were rationalized using a kinetic model, where 
charge transport occurs via a hopping mechanism between 
active sites delocalized over four base pairs. These results 
suggest that DNA—C61 are interesting candidates for the 
realization of nanowires for future molecular electronic 
applications. 


E EXPERIMENTAL METHODS 


Figure la shows an AFM topographic profile of a single DNA— 
fullerene (DNA—C61) molecule which was deposited on an 
a) c) 


39 nM On, 
4 


102 nucleotides 


Figure 1. (a) Contact-AFM image of a single DNA-fullerene 
molecule on top of ultraflat gold (cantilever spring constant: 0.4 N/ 
m; scan size: 0.25 X 0.25 um’). (b) Schematic representation of the 
DNA-—C61 molecule and (c) chemical structure of the C61—ssDNA 
binding group. 


ultraflat gold surface. As can be seen in Figures 1b and 1c, the 
DNA-—C61 molecules were obtained by the formation of an 
amide bond between carboxylic group of Cg, CHCOOH (C61) 
and the 3’-position of amine-modified single-stranded DNA 
(NH,-ssDNA) (102 nucleotides, 66.7% GC content). The 
specificity of the amide bond allows for the formation of only 
one C61 group attached to the 3’ end of one single-stranded 
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DNA.” The same procedure was repeated for the comple- 
mentary NH,-ssDNA before hybridization of both strands. 

After purification and characterization (Figures S1—S3), the 
double-stranded DNA-—C61 molecules were deposited on 
ultraflat gold substrates (Platypus Technologies, LLC; 
Madison, WI) and dried for 24 h inside a desiccator (see 
Supporting Information for details). The GC content of the 
DNA sequence was intentionally designed to favor a G-hoppin 
charge transfer mechanism, as suggested in the literature.” 
Moreover, fullerenes are known to form better electrical 
contacts with gold than e.g. thiols and amines,” hence 
allowing for C-AFM experiments at room temperature.” 

Our synthesis protocol is based on the work of Shin et al.,”” 
with minor changes. The compounds Cgg°CHCOOH (C61), 
N-hydroxysuccinimide (NHS), and 1-ethyl-3-(3-(dimethyl- 
amino)propyl)carbodiimide (EDC) were purchased from 
Sigma-Aldrich (St. Louis, MO). The DNA sequences 
(ssDNA) were purchased from Eurogentec (Liége, Belgium). 
The exact sequence is reported in Supporting Information 
Section 1. After functionalization and removal of the unreacted 
compounds through dialysis, FT-IR spectroscopy (Bruker 
Tensor 27, Pike Miracle ATR cell) was used to verify that an 
amide group was formed between the fullerene and the DNA. 
The next step was to anneal and hybridize the two ssDNA 
strands. In order to verify the occurrence of only a single band 
of modified double-stranded DNA-C61 without impurities, the 
samples were analyzed on a 2% agarose gel with a reference 
DNA ladder, the GeneRuler 100 bp DNA ladder (Thermo 
Scientific Inc.). These DNA—C61 samples (0.1 mM PBS) were 
used for C-AFM and STM experiments. In order to obtain 
single molecules of DNA—C61, the original stock was diluted 
in Milli-Q water, drop-casted onto a flat gold surface, and dried 
with N, flow, yielding ~10 molecules per pm’. 

This sample was left in a desiccator for 24 h before 
performing analysis by either atomic force microscopy (C- 
AFM, Bruker Multimode 8 AFM) or scanning tunneling 
microscopy (STM, Omicron Nano Technology, data shown in 
the Supporting Information). The C-AFM measurements were 
done inside a Faraday-shielded glovebox with protective N, 
atmosphere at room temperature and relative humidity below 
6%. The STM measurements were performed at room 
temperature under ultrahigh vacuum (base pressure below 5 
x 107"! hPa). 


© RESULTS AND DISCUSSION 


All AFM experiments were performed with a Pt—Ir AFM tip in 
a nitrogen environment (relative humidity <6%), using a 
Faraday-cage adaptation inside the AFM (Bruker Multimode 8 
AFM) to avoid electronic disturbances. A control experiment 
was done using C61, solubilized in Milli-Q water and deposited 
in gold (Figures S4 and SS). For each AFM measurement, the 
total height (2.6 + 0.1 nm) and length (39 + 1 nm) were 
verified to correspond to an individual double-stranded DNA— 
fullerene molecule (DNA-C61 of 102 bp) in its standard B- 
form and lying horizontally on the surface (see Figure la).”* 
Because of the drying procedure, one might assume that B- 
DNA can undergo a structural transition to A-DNA or to Z- 
DNA in the absence of water. However, this would result in a 
compressed or a supercoiled structure, which was not observed 
in our experiments.” This suggests that the low, but nonzero, 
humidity level in the AFM chamber is still sufficient to stabilize 
the B-conformation via hydrated counterions.”* 
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Figure 2. Schematic representation of the C-AFM experiments. (a) The Pt—Ir tip is brought on top of the DNA—C61 molecule at the starting 
position T1 (middle of the molecule) and voltages between —1 and +1 V are applied. I(V) curves are measured stepwise every 0.65 nm. (b) Change 
of I(V) curves when the AFM tip moves from T1 toward a fullerene end. The current cutoff at +5 nA is due to signal saturation. (c) I(V) curves 


obtained with five different molecules for three representative tip positions. 


C-AFM measurements were performed by positioning the 
AFM tip in the middle of individual DNA—fullerene molecules 
and applying a voltage between the tip and the gold surface 
(see Figure 2a). Voltage scans (— 1 to +1 V) were done at 30 
different positions on each molecule, starting from the middle 
and going toward one of the fullerene ends and finally out of 
the sample with a fixed translation step of 0.65 nm. Figure 2b 
shows the I(V) curves measured at selected positions along a 
given molecule. As seen in this figure, we observed a clear 
steepening of the I(V) slope when the tip was moved from its 
starting position (T1) in the middle of the molecule to the 
fullerene end (T30). The I(V) slope at T1 corresponds to a 
resistance of ~200 MQ and decreases to 20 MQ at position 
T30. We note that the apparent current cutoff at +5 nA is due 
to the limited current range of the C-AFM setup. The complete 
set of curves can be found in the Supporting Information 
(Figure SS). We also found that the type of base pairs in direct 
contact with the AFM tip seems to affect the I(V) slope. For 
regions rich in G-C base pairs, the I(V) curves are almost linear 
while they present a typical semiconductor aspect in regions 
containing mainly A—T pairs (Figures S6, S7 and Table $1). In 
addition, we also investigated the structural effect of desalting 
the DNA—C61 molecules, which resulted in very low currents 
compared to the standard DNA—C61 molecules (Figures S8, 
S9 and Table S2). 

Identical I(V) measurements were performed on 22 different 
DNA-—C61 molecules. Conductive behaviors similar to those 
shown in Figure 2b were observed in 19 out of 22 of these 
samples. To illustrate the relative reproducibility of our 
experiment, we show in Figure 2c the I(V) curves of five 
different molecules measured at the selected positions T1, T15, 
and T30. Although the measurements were taken over a period 
of several days and that there is a subnanometer uncertainty in 
hitting the exact starting position T1, the different measure- 
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ment present the same trend as the conductance increases 
strongly when the tip approaches the terminal fullerene group. 

The variation of the conductance observed in Figures 2b and 
2c with respect of the tip position clearly indicates that the 
charge transport occurs along the entire length of the DNA 
strand and not in perpendicular direction between the C-AFM 
tip and the surface through individual base pairs. In addition, 
experiments performed on freshly prepared samples using a 
scanning tunneling microscope confirmed the absence of any 
type of vertical tunneling mechanism through these molecules 
(Figures S10 and S11). The possibility of direct charge 
tunneling between the AFM tip and the surface was also 
tested using coherent charge transport simulations (see Section 
5 of the Supporting Information). These simulations show that 
the current intensities due to this mechanism would remain 
close to zero in the voltage range used here. This is due to the 
large band gap of the molecule and the very weak electronic 
couplings (0.1—1 meV) between the molecular orbitals of the 
base pairs and the electronic states of the electrodes. Thus, 
charges are more likely to propagate through the a-stacked base 
pairs (Figures S$12—S14), where electronic couplings are strong 
(50—150 meV), rather than to tunnel from the tip to the 
surface (Figures S$15S—S17). 

The large current intensities reported in Figure 2 confirm 
that the strategy employed here to deposit DNA—C61 on a 
gold surface preserves the B-DNA structure and does not 
disrupt the z-stacking between neighboring base pairs. An 
important point to note here is that as mentioned above, the 
shape of the I(V) curve changes when the nature of the base 
pairs that is in direct contact with the AFM tip changes. This 
shows that the high currents observed here are due to transport 
through the DNA z-stack and transport via salt atoms adsorbed 
on the DNA can be ruled out.”* 


DOI: 10.1021/acs.jpca.7b00348 
J. Phys. Chem. A 2017, 121, 1182—1188 


The Journal of Physical Chemistry A 


a) 


Fin (Yki 


Kur f 


ky 
~ Cas - A A, N ao zo 
e a O e 00. 
f 


e 


ksurf 


Figure 3. Representation of the kinetic charge-transport model. (a) The model involves multistep charge hopping between N active sites, each of 
them corresponding to several interacting base pairs. The injection and extraction rates Keo are here calculated following ref 38. The inter-base- 
pair rates k, and kç are calculated via a balance equation. (b) The voltage drop along the DNA is assumed to be linear. 


As an additional confirmation of the direct relation between 
the conformation of the DNA and its conductive properties, we 
have also studied desalinated DNA—C61 molecules. In these 
experiments we observed considerably reduced conduction 
values (Figures $8, S9 and Table S2). This is as expected since 
desalination is accompanied by a disruption of the B-DNA 
structure and a poor z-stacking along the strand.'> Hence, the 
presence of ions and transport along the z stack are indirectly 
related. 

To understand the origin of the large current intensities 
observed experimentally, we have modeled the charge transport 
along the DNA using the rate equation proposed by Lehmann 
et al. Figure 3 illustrates the assumptions and parameters 
entering this model (see Section 4 of the Supporting 
Information). 

As seen in Figure 3, charges move via a series of incoherent 
hopping steps through a sequence of neighboring active sites. 
An important assumption of the model is that in the absence of 
voltage drop all the sites present the same energy. This 
assumption, which might seem counterintuitive given the 
heteropolymeric sequence studied here, was made to provide 
the simplest possible picture of the charge transfer mechanism. 
Accounting for the exact sequence would require the use of a 
full quantum mechanical description of the system or the 
introduction of a large number of parameters in the model. The 
former is beyond the scope of this article, and we believe that 
the latter would only cloud the understanding of the charge 
transfer mechanism. The voltage profile along the molecule is a 
key ingredient of this model, as it defines the energy landscape 
for charge propagation. Following previous work,” we assume 
that the applied voltage drops linearly between the two termini 
of the molecule. The voltage profile is then governed by three 
parameters: @,;, and @,,,,¢ define the voltage drops at the tip and 
at the surface while a@,,,; is the total voltage drop along the 
molecule. The conservation of voltage imposes that Aup + Our 
+ mol = 1. Therefore, the voltage dependence of the current is 
determined by the interfacial charge transfer rates between the 
metallic contacts and the molecule kahat Following Polizzi et 
al? the charge injection rate from the tip to the molecule is 
given by 


A+ 5 + ayeV } 


kip(V) = f dx exp +: — T 
z 5 


kT x\-1 
«(= (1 +e”) 


(1) 
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In this equation, Cup represents the coupling strength at the tip— 
molecule interface and A is the reorganization energy 
corresponding to the charge injection process. The other 
injection rates have similar expressions, described in Section 4 
of the Supporting Information. The forward and backward 
charge transfer rates between different active sites along the 
molecule are determined by a balance equation 
k; = koe [mae V/C(N-1)ksT)] and k, = kget roe V/(2N-1)6T)) 


where ky is the charge-transfer rate between neighboring sites 
in the absence of applied voltage. These voltage-dependent 
transfer rates are used as parameters to solve a kinetic model 
given by the equations 


È = —(hip + kB + kP, + k apb tad (2) 
É, = — (kp + ke)P, + kPa + keP g1 
Py = — (kit + ke) Py ag kePy_4 + Kanee contact 


This system of equations is solved in the steady state, i.e., Py = 
0, and with the normalization condition P, + P, + .. + Py + 
Peontact = 1. Once the steady state populations of the sites are 
calculated, the currents flowing through the tip and the surface 
are obtained by 


I.(V) = = 0K pF ontack T kipP) 


Lwt( V) = —e(ksutPN T ko 


surf 


(3) 


ae) 


Note that the values of Iyp(V) and I,,,.(V) are identical in the 
steady state limit. In its present form, this model can only 
simulate two-terminal systems, and therefore we considered 
only the charge propagation between the AFM tip and the 
closest fullerene terminus. Hence, the calculated total current 
intensity is possibly underestimated. Because of the strong 
electronic coupling between the fullerene and the surface 
(Figure S15), its orbitals are hybridized with the electronic 
states of the electrodes, making the last base pair of the 
sequence the final active site in the hopping model. Thus, the 
interfacial rates k4“ represent the charge injection/extraction 
between the fullerene and the DNA. The electronic couplings 
supporting this mechanism are however weak due to the linker 
group binding the fullerene to the DNA backbone. 

Previous reports suggest that the charge density can be 
delocalized over up to six base pairs during its propagation 
along the DNA sequence.””** Therefore, we consider that each 
active site in the model corresponds to a short sequence of N 
neighboring base pairs. This delocalization reduces the static 
energy disorder of the active sites by averaging the energies of 
individual base pairs. 


DOI: 10.1021/acs.jpca.7b00348 
J. Phys. Chem. A 2017, 121, 1182—1188 


The Journal of Physical Chemistry A 


15 


L (nm) 


Figure 4. (a) Fit of I(V) curves for a given DNA—C61 molecule at three different electrode positions using the kinetic model. The calculated curves 
(solid lines) are in good agreement with the experimental data (shown as circles). (b) Conduction variation of this molecule as a function of distance 
L between the AFM tip and the terminal C61: Data are shown as black squares while the calculation is shown as dotted blue line. The gray square 
dots represent data for different DNA molecules at the same tip position along the double strand. 


Table 1. Values of the Parameters Entering the Kinetic Model Fitted to the Experimental Data Shown in Figure 4a 


ô (meV) (6) (meV) A (meV) fy (COM 
T1 65 25 85 10 
T15 65 25 85 10 
T30 65 25 85 10 


Cone (UO 6) ip Osarf ko (10"° s7!) Nyite 
10 0.1 0.2 23 13 
10 0.1 0.2 23 7 
10 0.1 0.2 23 1 


“The position T1 is model with 13 active sites but corresponds to 51 base pairs between the tip and the end of the molecule. This leads to about four 
base pairs per active site. Note that aside from the number of active sites, all the parameters are kept identical for the simulation of charge transport 


for different tip positions. 


However, the disorder of DNA molecules at room 
temperature leads to fluctuations of these energies. Hence, we 
assume that the energies of the active sites are distributed 
following a normal distribution centered at 6 (Figure 3b) with 
standard deviation (6). These fluctuations account for the slow 
modifications of the molecular structure. These modifications 
induce a change in the energy landscape that each electronic 
charge experiences during its transfer. The parameters entering 
the calculation of the current intensity (eq 3) were fitted to 
reproduce the experimental results displayed in Figure 2b. Since 
these experimental data were obtained for a single molecule, all 
model parameters were kept identical in the calculation with 
the exception of the tip position on the DNA—C61 complex. 
During the fitting procedure, the rate equation was solved M = 
250 times to average the effect of the static disorder. 
Characteristic fits are shown in Figure 4a, and the parameter 
values are summarized in Table 1. Our numerical simulations 
indicate that each active site comprises N = 4 neighboring base 
pairs, which is in good agreement with other recent results.”””** 
The delocalization occurs probably within the blocks of 
identical neighboring base pairs contained in the sequence. 
Note that these blocks do not contain exactly 4 base pairs, and 
the delocalization length should be seen considered as an 
average size of these blocks. 

As seen in Table 1, the fit to the experimental data leads to 
an energy value for the active site of 6 = 65 meV. This value is 
relatively low compared to recent studies that lead for example 
to 6 = 550 meV for G quadruplexes.** However, the lack of 
band gap in the experimental I(V) curves shown in Figure 2 
clearly point at a low value of 6. This value is a dictated by 
several parameters such as the nature of the electrodes and the 
structure of the molecule and therefore might vary significantly 
from one experiment to another. Our analysis also suggests that 
the static disorder is rather limited with a standard deviation 
(5) = 25 meV. For comparison, other studies***°-** reported 
variations in the energies of the individual base pairs of up to 
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100 meV. Note however that these previous estimates of the 
static disorder were obtained for a DNA hairpin in solution and 
not deposited on a rigid surface. This might explain the reduced 
static disorder observed in our experiment. The transfer rate 
between neighboring sites is here equal to ky = 2.3 X 10'! s. If 
we assume an inter-base-pair electronic coupling of Vp, ~ 50 
meV and a reorganization energy of Ap © 0.5 eV as given by 
previous studies,’”***° this value of the inter-base-pair transfer 
rate is compatible with an adiabatic charge transfer mechanism 
between the two delocalized sites (see Supporting Information 
Section 6). 

This relatively high value obtained for kọ gives an additional 
confirmation that the structural integrity of the B-DNA is 
preserved in our experiment. A disruption of the z-stacking 
between neighboring base pairs would result in lower the 
charge transfer rates. Finally, the transfer rates cy, and Cup are 
both equal to 10 x 10° s~!, and the reorganization energy 
associated with charge injection is 85 meV. These values, 
combined with simple calculations of the electronic coupling 
between the AFM tip and the DNA reported in Figures S16 
and S17, are consistent with a nonadiabatic charge transfer 
mechanism between the electrodes and the base pairs (see 
Supporting Information Section 6). 

We have finally studied the distance dependence of the zero- 
bias conductance G = dI/dVly_y obtained along the DNA—C61 
molecules. The experimental values of conductance were 
directly extracted from the slope of the I(V) curves measured 
for the 19 molecules and are reported in Figure 4b. The 
experimental data show a large spread but reveal a clear 
common trend as the conductance decreases with the distance 
between the tip and the terminal fullerene. Note that the 
deviations at distances exceeding 10 nm appear magnified due 
to the logarithmic scaling of the conductance axis. The values 
obtained for one specific molecule are highlighted by black 
squares. This shows the typical conductance variation observed 
during a single experiment. The theoretical curve of 
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conductance vs distance was then calculated with the multistep- 
hopping model presented above using the values of the 
parameters reported in Table 1. As seen in Figure 4b, there is a 
good agreement between theory and experiment as both show 
similar distance dependence. To confirm that charge transport 
occurs only along the DNA sequence, we have introduced a 
direct charge transport from each active site to the gold surface 
in our model (see Section 4 of the Supporting Information). As 
seen in Figures $13 and S14, such direct charge transfer leads to 
a much weaker distance dependence of the conductance and a 
pronounced nonlinearity of the I(V) curves. Both are in 
contrast with the experimental results, which suggest that direct 
injection plays only a negligible role. 


© CONCLUSIONS 


Our results indicate that DNA—C61 molecules are good 
candidates for the design of high-conductance nanowires. Our 
AFM characterization revealed that the structural integrity of 
the molecules was preserved during their deposition on the 
gold surface. The DNA then assumes its standard B-form with 
good a-stack interactions between neighboring base pairs. 
Conductive AFM measurements of the conductance of these 
complexes showed that DNA—C61 complexes are able to 
efficiently carry charges over large distances and with large 
current intensities. In addition, separate measurements on 
different molecules and on different days lead to similar results, 
illustrating the relative reproducibility of our results. These 
experimental findings were reinforced by a simple theoretical 
analysis based on a rate equation. These simulations indicate 
that an incoherent charge hopping mechanism could be 
responsible for the propagation of charges. The charges hop 
between identical active sites each containing four base pairs 
due to the charge delocalization. Despite its apparent simplicity 
compared to the complex sequence of the DNA molecule, this 
model was able to reproduce individual I(V) curves and also 
lead to a distant dependence of the conductance that is in good 
agreement with our experimental data. The results presented in 
this article show that DNA—C61 molecules provide an 
attractive alternative for the production of nanowires that 
could be used as connections in electronic devices at ambient 
temperature. 
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